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ABSTRACT 


The increasing weight and speed of carrier based aircraft are 
taxing the limit of conventional piston type energy absorbers which 
are used for arresting gear. Turbo-type absorbers have been proposed 
as an alternative. This study investigates a turbo-type energy 
absorber for an aircraft arresting gear under development by the 
Naval Air Engineering Center, A mathematical analysis, computer 
Simulation and performance Reeder on are given for each mode of 
operation, It was initially expected that the absorber could operate 
gm just two ae ee and reverse, but the analysis shows that 
- four distinct modes are possible, since not only shaft rotation but 
the path of fluid flow can be reversed. Theoretical performance 
predictions are also compared with test data from an existing smaller 
scale version of the Becks, The agreement is excellent. An 
epproximation of scaling effects on power absorption is also 
presented. It is concluded that the turbo-type absorber is basically 
adequate to meet the demands of carrier operation for the forseeable 


Devure. 
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I, INTRODUCTION 


Over the years, the speed, size, and weight of carrier based air- 
craft have steadily increased, thus demanding an ever larger energy 
absorption capacity from the shipboard arresting gear. Up to now, 
design improvements to meet these increasing demands have been limited 
to mere enlargement of the basic linear hydraulic energy absorber, 
which has been of the piston type. Modern aircraft, however, have 
joaded this system to its upper limit, and make any further development 
along these lines impractical due to weight and space limitations. | The 
need exists for a new concept which will minimize space and weight 
demands while providing the necessary increase in energy absorption 
Cepeacity. 

The solution proposed involves the use of a turbo-type rotary 
hydraulic energy absorber. In essence the fluid is put into motion by 
a suitably designed rotor and is Pisa joe ana pyotne 
resistance of the flow passage. The power input to the rotor is 
absorbed by the hydraulic resistance of the system, and is eventually 
dissipated as heat. 

The energy absorber is directly coupled to a cable storage drun. 
The drum profile is programmed to match drum rotational velocity to that 
required for optimum retarder performance. A device of this type was 
introduced in 1968, and is known as the M-21 energy absorber. ‘The 
initial design gives satisfactory performance, but is of rather limited 
capacity, thirty million ft.lbs. Therefore, a program has been under- 
taken to redesign the absorber to increase the capacity to fifty-five 


meiiion it.lbs. 
ee 








A development program was initiated by the Naval Air Ingineering 
Center (NAEC) (Ref. 1). A full scale test of the programmed drum 
coupled with the M-21 absorber was conducted to determine the 
feasibility of the overall concept and to obtain performance data 
needed in the redesign program. A schematic diagram of the proposed 
system is shown in Figure 1. 

A theoretical approximation of the flow conditions in a machine 
of this type was derived by Professor Michael H. Vavra (Ref. 2). To 
facilitate investigation of the effects of various design parameters 
NAEC formulated a computer program, From the computer analysis, a 
workable combination of design parameters was chosen such that torque 
and energy requirements were satisfied at the design point, the con- 
dition at which the absorber is expected to perform most often. 

Professor Vavra has also been commissioned to develop and test a 
flow model of the new design. This model will be useful for refining 
losses, checking for flow separation, and obtaining information that 
is not available from the mean streamline analysis. 

In this paper, a theoretical simulation and performance analysis 
of the energy absorber in all modes of operation is presented in texms 
of dimensionless parameters that facilitate comparisons of various 
designs and interpretation of model tests, The retraction problem 
is also studied, and preliminary estimates of the power needed to 
retract the programmed drum are given. 

A throttling system is incorporated into the design to aid in 
arresting aircraft of various sizes. The discussion and analysis 


presented here are concerned with the condition of open throttle only. 
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II, DESIGN VALUES AND CONSIDERATIONS 


The main features of the energy absorter design that is analyzed 
here are summarized in Figures 2 and 3, and Table I. 

The customary approach in the design of turbomachinery is first to 
optimize conditions with respect to a single specified design point. 
The resulting performance over a broad range of off-design conditions 
is then investigated to determine overall acceptability. The condition 
at which the absorber is expected to operate most often is naturally 
chosen as the design point. Specifically, a design torque of 
340,000 ft.1lbs. pete at a rotational speed of 76 rad./sec. 

-The maximum rotational speed the cable will impart to the drun, and 
eenscauently to the rotor, is about 85 rad.jsec. This iimit is fixed 
by the allowable landing speed of the aircraft, and by drum design. 

The mechanical power imparted to the rotor is continuously 
converted to thermal energy of the fluid. The temperature rise of the 
fluid is significant and imposes a constraint on the design. ‘The . 
temperature limit is the maximum allowable temperature the fluid may 
attain without causing cavitation. Obviously, heat must be removed to 
maintain equilibrium conditions. For this reason, provisions are made 
to bleed off a portion of the pace ene eealetd continuously and to 
replace it by fresh cool fluid. During a rapid series of arrests, the 
cooling problem is especially critical. 

The fluid used in the absorber Po eavet ay a elereod mix tiuee.. se lmEs 
fluid is acceptable because it remains usable over a wide range of 


temperatures. 
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Table I PARAMETERS OF THE A3SORBER 


14.61 inches 


Ra 

R, = 16.25 inches 
fe = 16.07 inches 
Ry = 14.93 inches 
b, - )2-0e inches 
be = 2,02 inches 
2a 4.18 inches 
rs = 2,02 arches 
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9 = 10.26 inches 


iiemper Of blades for rotor and stator: /2 
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Fluid properties: 
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A detailed discussion of the effects of various parameters and loss 
coefficients is undertaken in the analysis section. A few words, 
however, are in order here. The radius of curvature re of the inner 
contour of both the upper and lower passage must not be smaller than 
the blade height. This insures that the flow does not separate, and 
it also increases absorption capacity. Hstimates of losses in this 
machine are obtained through various empirical turbomachinery cor- 
relations as is shown in the analysis. Enthalpy losses in the passages 
are significant. Much of this loss is due to the bending of the 
passage with friction responsible for the remainder. However, a 


change in the coefficient of friction c, which expresses the direct 


its 
wall friction losses in the passages has little influence on overall 
performance. A value of Co = 0.016 was chcsen since it gives best 
agreement with M-21 test data. 

The blade shapes first optimized by Professor Vavra from a 
performance standpoint were found by NAEC to be unacceptable on a 
stress basis. Blade thickness, spacing and angles were accordingly 
revised by NAEC, and a blade of sufficient strength was found; its 
characteristics are listed in Table l. | | 

The reliability of this analysis can only be judged experimentally 
by correlation with test data. The theory correlates very well with 
actual test results on the M-21 design. From this it appears that the 
theoretical analysis is reasonable and that the new design should be 
workable, if not optimun. 

Beer caevies operations often require a series of rapid and 


repeated arrests. The drum must be retracted thus reversing the 
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rotational direction of the rotor. This process is investigated to 
obtain preliminary estimates of the power required to retract the 


absorber. 
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Jil. ANALYSIS OF THE PROBLEM 


A. MODES OF OPERATION 
The energy absorber is capable of operation in four modes. These 


modes are summarized in Table II. 


Table It MODES OF OPERATION 


Flow Direction shaft Rotation 
Normal Reverse 
Normal Mode 1 Mode 2 
Reverse Mode 3 Mode 4 


Mode 1 corresponds to design point conditions. In this mode, the 


A 1 a: 
J uv 


y  2igeolee a GliGs vac 


fluid passes successively tnrougn stations Shi; 
in that order. Note that the flow leaving the LOLOL al Sitaaeton 2 ys 
moving radially outward. See Figures 2 and 3. For fluid circulation 
in this sense, conditions at stations 2 and 4 are particularly sig- 
nificant. At the rotor exit, station 2, the direction of the flow 
feu respect to the rotor, is assumed to remain constant at all flow 
rates. At stator exit, station 4, the direction of the flow with 
respect to the stator is likewise assumed to agree constant. While 
this assumption of fixed relative and absolute flow directions is not 
exact, it represents an excellent approximation when the blades are 
closely spaced as in the present design. 

| In Mode 2, the sense of the fluid circulation is exactly as 


described above. However, the direction of shaft rotation is reversed. 
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Naturally, stations 2 and 4 are still governing in the sense that flow 
directions at these stations are known. 

In Mode 3, the shaft rotation is the same as for Mode 1, but the 
direction of fluid circulation is reversed. Consequently, the fluid 
passes successively through stations 2, 1, 4, 3 and back to 1, in 
that order. Note that the rotor exit is now at station 1, and the 
direction of flow at this station is radially inward. In this case, 
stations 1 and 3, are controlling. The relative flow direction at 
rotor exit, station 1, and the absolute flow direction at stator exit, 
station 3, are now assumed to remain fixed. 

In Mode 4, the sense of fluid circulation is the same as in 
Mode 3, but the direction of shaft rotation is reversed. Of course, 
stations 1 and 3 still govern flow direction as in the previous case. 

The same set of eauations describe Modes 1 and 2, but a separate 
set of equations is needed to describe the reverse flow Modes, 4 and 4. 
Note that only two sets of equations are needed for the four Modes, 
and each set is governed by the flow direction relative to th blade 
eXits, Two separate computer programs are needed since one program 
will be used to solve each set of equations. Program Normal finds the 
solution for Modes 1 and 2, and <n Reverse finds the solution Poms 
Modes 3 and 4. Both programs will be discussed in the appendix. All 
quantities and results discussed throughout the paper are obtained 
from these programs. 

The velocity diagrams for each mode are presented in Fisures 4am), 
6, and 7. All angles and velocities are drawn to scale. 

The thermodynamic processes in the absorber ae shown on the 


enthalpy-entropy diagrams, Figures 8 and 9, The process for nornal 
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Figure 4 Velocity Diagrams for Mode 1 Operation at 
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Figure 5 Velocity Diagrams for Mode 2 Operation at 
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Scale: O.5inch=100 ft/sec 





Figure 6 Velocity Diagrams for Mode 3 Operation at 
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Figure 7 Velocity Diagrams for Mode 4 Operation at 
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circulation, Modes 1] or 2, is shown ome se whl eoat or 
reverse circulation, Modes 3 or 4, is shown in Figure 9. The 
direction of shaft rotation affects the exact quantitative enthalpy 
values. The graphs, however, are merely qualitative. 

A "total head" H (per unit mass) is used which differs from the 
standard total enthalpy ne This is discussed in greater depth in the 
analysis section, but a few words are needed at this point to clarify 
the discussion. 

The fluid is considered incompressible which allows the thermal 
and mechanical effects to be uncoupled. A "total head" is 


conveniently defined as 


P 2 
p p 2 


Henceforth H is referred to simply as head. The head and total 


‘enthalpy are linked by the simple relation 


H = H, -u 


where u represents internal energy per unit mass. 


B. DIMENSIONAL ANALYSIS 
Results are presented in terms of appropriate dimensionless 
quantities. The use of non-dimensional terms simplifies the cor- 
relation of data, the comparison of various designs, and the inter- 
pretation of model tests. 
' The three quantities chosen as consistent dimensional reference 
parameters are the radius of the rotor outlet Ros the fluid density P, 


and the shaft rotational speed W or W » whichever is appropriate. 
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ay 








Using these parameters, one defines the effective Reynolds number 


for this design to be 





where v is the kinematic viscosity. And the dimensionless total head is 





In the discussion of results, it is convenient to define the flow 


coefficient ¢ in the following specific way, namely, 


» - Ll 


WR, 


where the algebraic signs of w and ¢ are always positive for normal 
shaft rotation and negative for reversed shaft rotation. This sign 
convention for shaft rotation is maintained regardless of the direction 
Sere fluid circulation, hence the use of the absolute magnitude eral 
in the definition. 

In the actual computer programs NORMAL and REVERSE, however, some 
what different sign conventions are more convenient. All equations are 
vei FEN in such a way that the three quantities Viyo WE, and ¢ are always 
positive. For this reason the four modes are computed as separate cases. 

For comparison of various modes and designs the dimensionless torque, 
power and rotational speed coefficients are respectively 
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wax Bove 


30 











lp = 

ioe By 
qh ia = 

MAX 


The detailed analysis considered later also uses various 
dimensionless quantities for convenience in the derivations. The 
parameters used to facilitate the non-dimensionalization in that case 
are not necessarily the same as introduced in the present section, and 


the separate cases should not be confused. 


CG. ENERGY BALANCE. AT EQUILIBRIUM 
In general the rotor adds energy to the fluid. In traveling 
“through the passages and returning to its initial state the fluid 
dissipates this mechanical energy through various losses. Under 
equilibrium these losses exactly cancel the original energy input. 
it turns out that the energy input of the rotor can be expressed 


mimcerms of the flow coefficient ® in the form 


= g 
bh Hoy Ao 


where A and B are true constants. The losses can be expressed as 
y SAC enn aac ce 
AHL oss a 


where a 3 and C! may be considered constant to a first approximation. 
Actually these quantities are weak functions of $, a fact which must 

be considered in the exact solution. The energy relations which fix 

the equilibrium state are shown in Figures 10 and 11. Figure 10 


represents the two modes which correspond to normal flow direction. 
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Figure 11 shows the two modes for reverse flow. These results are 
drawn to scale, except that the parabolic lines are only approximate. 
The final equilibrium points, where the losses balance the energy 
input, however, and the relative distribution of losses at these 
points are in correct scale, 

Initially, it was thought that only two of the above four modes 
would prove to be physically possible, one corresponding to normal 
Shait rotation, the other corresponding to retraction. Surprisingly 
enough, the present analysis has demonstrated that operation in all 
four of the modes is theoretically possible, This point is clearly 


tiemeated in Figures 10 and 11. 


D. MODEL TESTING 

Professor Vavra has proposed that a model testing program be under 
taken. The model will not be a true dynamic model, but will be a purely 
static device that uses air as the working fluid. Construction of the 
model will constrain the air flow to the normal direction only. However, 
the scale model should be useful in checking for flow separation 
through the passages and refining the loss coefficients for the 
analysis of the absorber. The effects of the tangential component of 
velocity will also be investigated. 

The discussion of the above static model is not within the scope 
of this thesis. However, a hypothetical model is considered which is 
a true dynamic model of the actual full scale device. The size of the 
hypothetical model is taken as identical to that of the static model. 
The model is scaled to 0.569 of the size of the actual absorber. In 


both cases the fluid is air. When expressed in dimensionless terms, 
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the performance of the hypothetical model is basically the same as that 
of the full scale machine, except for Reynolds number effects, The 
full scale machine operates at a maximum Reynolds number, as defined 
above, of 4.88 x 10° while the maximum Reynolds number of the model is 
Be ox 10°. The results for this model are presented along with the 


results for the real absorber. 


BE. ASSUMPTIONS AND LIMITATIONS OF THE ANALYSIS 

A one-dimensional mean streamline analysis is obviously only a 
crude approximation of the actual performance of any turbomachine, 
but is a starting point necessary in most development programs. ‘The 
validity of the assumptions used are very important to the reliability 

of the analysis. The primary assumptions used in this report are 
ettriced below. 

The fluid is treated as incompressible. This assumption is 
excellent as long as the provisions to maintian the fluid at a 
temperature lower than ae cavitation temperature are adequate. This 
conveniently uncouples mechanical and thermal effects and allows a 
form of “total head" H to be defined which is independent of 
internal energy. 

iiecne velocity distribution across every section were turly 
uniform, the meridional component V at. the mean streamline could be 


il 
computed simply from a continuity relation of the form 


where Q is the volumetric flow rate and Any is the meridional cross- 


section area. The actual velocity distribution, however, is non- 


DD 








uniform because of the effect of streamline curvature, boundary layers 
and so on. To allow for these effects, approximately, the actual 
meridional velocity Vu 1s assumed to be governed by a relation of the 


form 
Nu = %3 i 


where. the so-called blockage factor op is treated as a constant. The 
blockage factors were extrapolated from M-21 test data. Their reliability 
in the reverse flow modes is doubtful, however, since all test data 

were obtained only from Mode 1 operation. An extension of these data 

to all modes would be very useful. 

The present analysis treats operation of the absorber as a steady 
state phenomenon, Strictly speaking, the operation is actually 
unsteady. However, the arrest and retraction processes require a 
Erericient number of drum revolutions that the steady state assumption 
pmevtaes a fair approximation. 

A major assumption is that the relative flow direction at the rotor 
outlet and the absolute flow direction at the stator outlet are 
Sem-vant. Hor the case of normal flow this is an excellent ap- 
proximation. However, the validity is reduced for the reverse flow 
modes. in Modes 4 and 4, the flow encounters a sharp leading edge, 
and exits past a blunt trailing edge. The implication is that the 
flow is not as well guided as in the normal flow case. Thus the 
actual turning angle of the flow is undoubtedly less than predicted 
by the theoretical analysis. The power absorbed in both reverse flow 


modes would consequently be less than predicted. 
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Losses, in general, are approximated by various empirical cor- 
relations and are reasonable as shown by test results. They are 
discussed in depth in the detailed analysis section. 

Incidence losses, however, pose a problem. The incidence loss 
approximations introduced in the analysis section are reasonable as 
long as the incidence angles remain small. This is obviously not 
the case for the retraction modes as shown in Figures 5 and 7. The 
approximations break down physically at these severe incidence angles. 

It was thought earlier that two of the four theoretical modes 
might show themselves to be physically impossible. For example, a 
large incidence loss might conceivably preclude the attainment of the 
necessary energy balance. This has not proved to be the case, however. 
It appears from the present analysis that operation in all four modes 
remains possible physically. It is true, of course, that the predicted 
performance for the two retraction modes, Modes 2 and 4, will be 
inaccurate because of these large incidence effects, but this in itself 
does not eliminate the existance of these modes. 

At extreme angles of incidence, a distinct possibility of stalling 
exists. In normal operation, this limits the attainable power absorption. 
In a retraction mode, however, stalling would not be totally undesirable. 
If the rotor stalls, the energy absorbed by the device is greatly 
reduced, and the power necessary to retract the rotor is minimal. 

All of the foregoing assumtions are thought to be reasonable. 
Refinements of various points would be useful, especially 10m lange 


incidence angles and blunt trailing edges. 








IV. RESULTS AND DISCUSSION 


The discussion of the problem up to this point has already included 
a number of significant results reached in the analysis as shown in 
Figures 4 through 11. 

Additional detailed results for all modes of the absorber are 
presented in Figures 12, 14, 14, and in Table III. The results for 
Mode 1 of the dynamic model are also included. 

The main results of this investigation are summarized concisely in 
Figure 12. This is a logarithmic plot of the dimensionless power Il 


as a function of the dimensionless shaft speed w/e, Performance 


TAX’ 
-of the proposed absorber, of the hypothetical model, and of the 
existing M-21 machine are compared on this graph. Note that the 
theoretically predicted performance of the M-21, version as shown by 
the continuous line is io excellent agreement with the actual measured 
performance as shown by the data points. 

The actual horsepower is plotted as a function of absolute shaft 
speed for all four modes of the proposed absorber in Figure 13. ‘The 
power absorbed during arrest and the power expended during retraction 
are presented for a range of possible operating speeds. 

The maximum dimensionless power is ae shown as a function of 
Reynolds number Ra im Pisureria. Uaws graph, although based on only 
three computed eee does correctly portray, approximately, how 


scaling affects power absorption. It would be a simple matter to 


calculate additional points on this curve from program NORMAL. 
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As mentioned earlier, the actual existence of four stable modes is 
a surprising result. The original expectation, based largely on 
intuition rather than analysis, was that only one mode would be possible 
for each direction of shaft rotation. But the absorber is an unusual 
machine which tends to confound intuitions based on experience with more 
conventional devices. Consider for example the fact that losses are 
necessary for the machine to absorb the mechanical energy input. The 
smaller the loss coefficients, the higher the final equilibrium flow 
rate will be, and the greater the energy absorption. Thus we are led 
to the seeming paradox that the lower the loss coefficients the higher 
the energy lost! 

The power curves of Figures 12 and 13 are interesting and 


Significant in that all curves have nearly constant and nearly equal 


cr 


slopes. The power can be expressed in the form 


iq 


a Y tax a) 





where the slope of the curve is roughly 17\=3 for all cases. 

The horsepower curves of Figure 13 indicate that Mode ig bhe 
design point mode, absorbs the greatest amount of energy. Mode 3, 
however, is a very close second. Refinement of the analysis should 
lower the predicted absorption of Mode 3 somewhat. Nevertheless, the 
favorable comparison of this mode to the design mode is interesting. 
Figure 14 also shows that the power needed to retract the rotor is 
reasonably small for Modes é and 4, especialiy in comparison to the 
power absorbed during arrest. 

The Ese of which mode the absorber will actually operate in 


under various conditions of use is interesting, but by no means 
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completely clear at this time. It is thought possibly that if the 
absorber be started from rest with shaft rotation in the normal 
direction, Mode 1 siamanccout: If retraction follows immediately, 

the retraction mode will probably be Mode 2, This is based on the 

fact that during arrest the flow will build up a large momentum which 
will tend to persist in the same sense during the subsequent retraction. 
if retraction be started after the fluid motion ceased, however, the 
resulting mode of operation remains uncertain. 

Some method of sensing and/or controlling flow direction would be 
useful to insure that Mode T was occurring on all arrests. Mode 3 
operation might not have enough absorption capacity to handle large 
eirerait. 

The general question of which mode occurs is an area for further 
study. Testing of the actuai absorber with the idea of definitely 
determining the operating modes under various conditions would be 


extremely useful. 
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Ve. CONCLUSIONS AND RECOMMENDATIONS 


The development of the problem and the exact solution of the 
two sets of equations describing normal and reverse flow conditions 
clearly show that four stable operating modes exist. The analysis 
adequately predicts the absorption capacity of the proposed device 
for Mode 1. However, the analysis is not as accurate for the 
reverse flow or retraction modes. 

It is recommended that the retraction problem be further studied 
to refine loss coefficients and to more accurately predict the power 
necessary to retract the absorber. It is also recommended that 
further study be made into the question of predicting which monies 
occurs under various operating conditions. This study would 
encompass the unsteady phase of the problem. Actual testing of the 
absorber should also be conducted with the intention of determining 
which mode will occur. A method of flow control or sensing would 
also be useful to guarantee proper performance in that Mode 4 
absorption capacity may be insufficient in certain cases. 

This turbo-type of energy absorber appears to be completely 
adequate to meet the demand of carrier operations now and in the 
forseeable future. The potential for future development and 


enlargement is ample. 
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APPENDIX A 


DETAILED ANALYSIS OF THE ENERGY. ABSORBER WITH THE 
FLUID CIRCULATING IN THE NORMAL DIRECTION, 
MODES 1 AND 2. 


A one dimensional mean streamline analysis after Vavra (Ref. 2) 
will be developed with the following assumptions: 


B 


Constant 


I 


2 


, Smeons vans 


Uniform flow at stations 1 to 4 except for blockage factors. 
Blockage factors are independent of the flow. 
The fluid is incompressible. 


oteady state conditions exist. 


1.0 Torque and Velocity Relations 


The torque T, exerted on the fluid by the rotor is positive 


R 


when acting in the direction of rotation. 


T =o (RV, - 8, Vy] — (1) 


where hm is the mass flow rate. 


o A Jie 


p(2n R, be) Yuu, 55, (2) 


m 


I 


I 


o(2n Ry bp) ‘Mo Ee 


The quantity 27m R b, is the flow area, V and Vo are the meridional 


R m1 2 
velocities at the mean streamline, and Koy and Kao are blockage factors 


which are assumed to be independent of the flow. 
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By introducing a flow coefficient 


5 = M2 _ le (3) 
: Us, WR, 


all necessary velocity and thrust equatiors can be written in terms of 
known geometric parameters and the flow coefficient. An iteration 
technique is later employed in program NORMAL to calculate % The 
following velocity relations are developed now in order to simplify 
the energy relationships which follow. 

The sign convention used is that all peripheral velocity com- 
ponents are positive when in the direction of rotation. From Figures 


4 or 5 simple geometric relations are obtained which become, using 


eq. 3. 
W, = Vyp/cos B, = OR, 6/cos Bs (4) 
Woo = Vy> Tan B, = OR, $ Tan B, (5) 
Vio = Wi + OR, = wR, [1 + $ Tan B,] (6) 
Also 
Vio = Nhe Tan a, = OR, 6 Tan ,, ae 


An expression for o, in terms of the known relative exit angle can be 


obtained by equating the previous two equations, or 


Tan a, = Tan 8, + 1/8 (8) 


similarly at the rotor inlet 


Wo = Vy, Tan B, (9) 
= : 0 
ea (10) 


4? 








The mass flow equation and the definition of the flow 


coefficient yield 


on aro 


2 
Vv. =tor—- — (11) 
M1 2 R, Kay 
adic. 
R 
ae 
lee = $ WR, a —— Tan PR, + OR. (12) 


1 1 


The mass flow equations must be extended to include flow 


trouch the stator. 


m = 21 Ry b, Viyo Kao OTT Rz b. Vu Kes 


an R, bs P Vy Kpg 


(13) 


From the above relations, the meridional velocity equations are developed. 


“wr _ 72 *pe (14) 
Vo By Bay | 

“ya _ "2 po Pr 

Vino a Kz ie (15) 
Yyq _ "2 Spe Pr | 

Yo By Key 5 . (16) 


1.1 Frictional Losses in Upper Passage 


Let Mey be the frictional moment from forces along the walls 
of the upper passage from station 2 to 3, acting on the flow opposite 


to the direction of rotation. 


mR, Viz = mR, Vip - Mey (17) 
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setting Me = ie mR, Nei 


and ee (1 - r) 


eq. 1/{ becomes 


Ato 
%u3 ~ "ht By Yue 
From the mass flow equation 
Ven = GR, 2 Be OR 
M3 2 oe Kez be 


7.2 Prictional Losses in Lower Passage 


(18) 


(19) 


(20) 


(21) 


With M 4 being the frictional moment caused by forces along 


f 


the lower passage walls from station 4 to 1 acting on the flow opposite 


to the direction of rotation, the eauations for the lower passage 


are developed in a manner similar to the previous development. 


Meg = 


"hee 
“y 
Vo = "he Neal R, 


R dD 
eT 
me 4 Kea S 


t 
ab 
! 
> 


From Figure 4 it can be seen that 


Wan oO 


Von = "Mg 4 


Using the previous three equations gives 


R K. bd 
2 2 R 
Ne = T] é0R— ——- — Ten a 
WI ML 2 iy Bea be 4 
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(22) 
(25) 


(24) 


(5) 


(26) 


(27) 








From Figures 4 or 5, one can relate the absolute inlet angie 


to absolute velocity components as 


Tan @, = mel (28) 


Using the relations developed for the inlet velocity eqs. 11 and 27, 


the absolute inlet angle is 


lee 
Tan a, = BP ven a, (29) 
17 we % a 
Similarly the relative inlet angle is 
Veer rt 
gan p, = St 
M1 
WR 
= Tan = = 
1 Vo 
Subsvituting for Va yields 
Dee Ke eee 
R 1 1 
ines, = 1, = = Tana, = rae (30) 
1 "he BS Key %) Kao 


The relative and absolute velocities can now be obtained from 
their respective components by employing the following general equations: 


2 2 2 


V = Vy ave (aay 
Z 2 2 
Wav, +W, (32) 


ee 


= 
| 
<j 
we 
© 
O 
wn 
TW 
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1.5 Torque Relations 
The torque as defined in eq. 1 can now be expressed in 
terms of known geometric parameters and the flow coefficient. Using 


Paes. 51, 0, andui2 anecaq.. 1 ylolos 


b 


5 R my 
T, = an R, © Xoo @ Lt {fan 8o-Th,, = r= cE nan 045 — 


The non-dimensional torque coefficient 1, is now introduced as 


a 


b 2 
an(R,,)° “Ep Kyo 


S 


i= 


DL 


iv nic {Tan B, - ty = 5 Tan as (35) 


1 aA} alee 


The torque qs which the stator exerts on the fluid is positive 


wen acting in the direction opposite to rotation 


mT =m @, V,, - 8, Vy) (36) 


substituting the values for the whirl components gives 


iD 


: ee 
5OR\ 2 Pare a8? 
= = emR,, z) p Ww Kellie $(1+9TanB,,) - ¢ eo) Ke Tan 2h 


Again introducing a non-dimensional torque coefficient 


T 
S 4 4 
Sak 2 
em ae eeu? 
2 \B) 2 Nop fa 
pa 
2 R “pe 
i MoD = G os a, | 


» 





The general torque equation is 


tS a re 7 Ming ~ a 


The torque exerted on the rotor by the fluid minus frictional moment 
losses 1s equal to torque exerted on the fluid by the stator. 
The dimensionless torque coefficient introduced earlier is 


again given as 


T 
oe Wie 2 


2.0 HKnergy Relations 


The torque and velocity equations developed in the preceeding 
section remain in terms of the unknown flow coefficient $. In this 
section an algebraic energy equation will be formulated and solved 
to yield ¢, 

for the sake of this development, a total head term H will be 


defined which differs from the standard definition of total enthalpy 


Hn which is 
2 
7 . 
le te Se (39) 
Or 
2 P 
V ae 
Ete eae = Wat op (40) 


where h is the static enthalpy 
p is the static pressure 
uis the internal energy 


and Pa is the total pressure. 
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Since the fluid is treated as incompressible, one can uncouple 
the thermal effects from the mechanical effects, and define a 'total 


head! H which does not include thermal effects. This new term is 


denoted by 
ie z 
ta pv 
eS == = = tea 1 
gees (41) 


Obviously the enthalpies are related by 


(Hy -u) =H (42) 


This new definition of total head H is useful in that constant H and 
constant Pa lines are coincident. Throughout the remainder of the 
analysis, useful total head H will be referred to as simply head. Note 
- that this total head has the units of enthalpy (per unit mass). 

The process undergone in the absorber is now represented in 
Figure 8. Obviously the process does not returm to the original total 
enthalpy Hn? a ecrnee does return to the original total pressure 
and head. 

From Figure 8 the energy rise across the rotor is shown to be 


equal to the energy losses across the stator and through the passages, 


or 
AH, - AH - AH, - AH, = 0. (43) 
Therefore, 


(Hy - Hy) - (Hy - HZ) - GH, - Hy) - (Hy - Hy) = 0 (44) 


where all the above terms are written as positive quantities. 
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The four terms of eq. 44 will be derived in terms of the flow 
coefficient and the known geometric parameters as in Section 1.0. Kq. 
4h is then straightforwardly solved to yield $. Program NORMAL is 
employed to accomplish this task. Since various loss parameters 
depend on ¢, an iteration technique must be used. 

2.1 Conditions in Rotor 
The incidence angle of the rotor blading i. is defined as 


R 


or merely the angle between the tangent to the mean camber line at the 
leading edge and the relative inlet velocity. The incidence angle in 
is seen to become more negative in the direction corresponding to a 
greater lift coefficient for the blade. 

The relative velocity W) can be resolved inte tyo components 
respectively tangential and normal to the mean camber line. The energy 
associated with the tangential component Ww, is assumed to be largely 
recoverable. The energy apecaead with the normal component is 


assumed to be totally lost. 


The useful relative velocity at the rotor inlet is then 


W, =W, cos in : (46) 


and the velocity component lost due to incidence is similarly 


Wis = Wy San Lp (47) 


The basic energy equation of turbo-machinery is 


Hoo ~ Hoy = Yo Vo - UL Vy (48) 
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which upon expansion becomes 


V ey eile V at aioe 
z Se 2 2. _ Al i: 1 (49) 
Hn Hoy 7 2 9D 
Using static enthalpies, regrouping, and cancelling terms in eq. 48 
gives 
Wo - Up WS - he 
+ ————_—— = rr 
h, - h, + 5 (50) 


Substituting for the static enthalpy and rearranging yields 
P ee p ee U ef U 
p 7 0 2 2 2 1 
Upon resolving W,°/2 into its components and regrouping further, one 


obtains 








- (A+ +24 )- (uy - a - =) (2) 


This leads naturally to the definition of the “equivalent relative 


total pressure” as 





2 
Pe p W 
Be _ (= + = 
5 7 - ) (53) 
and 
ne a 2 
*El = = + er aE a (54) 
p p @ 2 


and to the definition of the rotor blading losses as 


ve iE - f 
ae oe (55) 


AP - 
= [M2 


ot 








It also proves useful to express this loss in terms of a dimensionless 


coefficient defined as 
- P 
ial E2 
elt) (56) 
LW 
a 2 
Recalling that 
H =H-u_ (42) 


and using the basic energy equation, eq. 48, one obtains the energy 


fee across the rotor in the form 
_ Pl 2 } 
aS - H,) = Up Vio - UL Va - {sin ° aoe oom (oe) 


This equation can be conveniently non-dimensionalized by 
dividing bv nee The velocity exoressions developed in Section 1.0 can 
not be utilized to simplify eq. 57, and the dimensionless energy rise 


across the rotor (IIA H),, can be written as 


7) H, - Hy (sin® oy) R, \2 
fo, eh 1 - (BY 
AH‘ R oR 2 2 Ry 


+ § {Tan Bp Es BS Tw Tan a, ) cosa} (58) 


re 


in-i R.2 2 bo 2 
#454 @ [E+ Ewma)] 


1 a 
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where the dimensionless head coefficient is defined as 


(ork 
Ut = oa (59) 
w R 
e 
The loss coefficient Yp of eqs. 56 and 58 is a function of the 


deflection 48, where 


aE (60) 


This will be further discussed in Section 3.0. 
2.c Performance of Rotor Blading 
The rotor is designed for the maximum useful enthalpy rise 
Necross the blades. If this were a frictionless process, the useful 
energy rise would naturally be greater than in an actual process 
involving friction. Therefore, the rotor efficiency is defined as 
He = 
Ta ; i (62) 
2 Th 1 
where the total heads H, > Hy» and Aor are shown on the H-s diagram, 
Figure 8. Of course the pressure rise through the rotor is related 


very simply to the corresponding energy rise according to the 


relation 
P 


He — (62) 
0 


D9 


The numerator of eq. 61 is obtained from either eqs. 57 or 
58. The denominator of eq. 61 is similarly obtained, but for the 
theoretical case both the incidence and blading losses are non-existent. 


Tt follows from eq. Somumee 


(i Ning fom 7 1) 
Atm/R =~ 2 _ 2 
Ww R, 
ti (63) 
eel ca® {Tan Bs - ty F 5 Kah Tan ay 
is the theoretical dimensionless energy rise across the rotor. The 
efficiency can obviously be written as 
m1 (64) 
IT AH va, r 


(Pama (65) 


————— (66) 


Recalling that 


| eee (3) 





m= emp R, bo, Var Kee (2) 
and 
Aer 
eS (59) 
W Ro 
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the dimensionless power and head coefficients can now be related as 


p= fon(G8) 


The power is related to the torque by the simple expression 





(67) 





19 Te (=) (68) 


2.3 Energy Losses in the Upper Passage 
The energy loss in the upper passage, assumed as only a total 


pressure loss due to friction between the rotor and stator, is expressed 


as 
Au, =H, -H,2-2 - 2 (69) 
Or - 5 
oP eb)- Ged) (0) 


idee Irteuuenless precess 


Va 
Hy = (Hyg y = a tS (71) 
and from eq. 1/ with the frictional moment Ma, equcimavOEment 
5 V3. 
(Vig)eq = Yuo R772 (72) 








2 (73) 


61 








Similarly for an actual process 


2 2 
D V V 
eS M3 u3 
Hp = ri a lee (74) 





the energy loss can now be written using the two previous equations as 


Pony ~ P2 al ‘ 
H, - H, = +H —3 + ( -1) 3 


2 ~ Hy ; oe (75) 


a 


"yh 


A convenient dimensionless pressure loss coefficient can now 


be denoted as 


Pp pe 
¢ = = | (76) 


(CB) 


r 
where C wilt be shown to depend primarily on the radius ratio 108 


the upper passage in Figure 2. This is similar to energy loss 
‘coefficients developed for bends in smooth pipes. 


Eq. 75 is non-dimensionalized, and using eq. 76 one obtains 





2 2 
Cc UV. Vv 
Set M3 al iL io u3 
(WAH) = 3 2, 2 ae (= 1) 2, 2 (77) 
re Yo et 


Using previously developed expressions for Vr one obtains an equation 


similar in form to eq. 58 or 


Ut 
Oo 


(1 AH), ona 


lL 


2 
GB ogrmet 


~~ 


Lp Dawe jc =) 5 

& he , 

: : - e z) on (2) ao Bo} 
3s ong e Ro 


4+ 
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2.4 Performance of Stator Blading 


From Figure 4 the stator incidence angle i. is 


i = @-—-ada_bB 
“3 “55 (79) 
t 
and the useful inlet velocity Vz is 
t 
Vz = Vz cos i, (80) 


The total head ahead of stator blades is given by 


Z 


De Ve Pa 
WM ote Sa 81 
b p Z p (0) 
while the useful head is 
fine t 
Ce eee 
' 5 3 eS 
5 p 2 p Kee) 


The above three equations combine to yield 


| (Vz')" = (v,)° 
a, = He + 5 
i (83) 
= H, - a3 (sin-i ) 
3 2 Ss 


similar to rotor blading losses, stator blading losses will be defined 


as 
! t 
yo 3d (84) 
Ss nee 2 
5V, . W,*/2) 
or 
2 a 2 
V V V 
_ ae ee ae (5) 
, = Ae os ao  = tke (sin i) + - oe o ; 








subtracting H, from both sides yields 


1 
2 eC 


wg aes (255,02 ee — 
te ee (86) 


The same non-dimensionalizing technique used in Sections 2.1, and 2.3 
is again employed and the energy loss in the lower passage is now 


available in the general form as 


2 2 
H, -H Isl= cael R 
ee Se) antsy 
w R w R © 


ef ey 


+ 4 | (sin® il 5) 2 Tarn 3, | 4 3° {(sin® i ) ao 8 a + Tan* 8 s 
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The above equation is actually 


es a Fy - H 
1 Ss ges ec 
pe 2 Zane Zee. 


W Ry Ww Ry, Ww Ry 


4 1 2 1 2 j tl (88) 


‘See el eee — a — See) ee eee —= C So) Gea 
e 


which gives the energy losses in the lower passage in terms of the 
energy changes throughout the remainder of the system. The next 


step is to find the energy loss in the lower passage and solve eq. 


44. 


64 








Again the loss coefficient Y, is 2 -funcrilon ol derlectien and 


the aspect ration, where 


AB = 4, - a (89) 


and 


en (90) 


2.5 Losses in Lower Passage 


The development of the energy loss in the lower passage 
closely follows the technique of Section 2.3. 


The energy loss in the lower passage is 


p v,° p v, 
Z epee ee peel eel 
A, =H, -H, => +3 cao) (91) 
for a frictionless process 
R 
a = 
(Viwdon =" B (92) 


1 


or relating theoretical to actual velocities yields 


<4 


(V4 en = Th, (93) 
Therefore, 
Ze 
Hy = (Hy) = sh + sh +3 ae) 94) 
while 
H, = - i ha " Yt , (95) 
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the pressure loss coefficient is 


2 
(D1 )a, 7 P4 ‘oe 
ame = Cy = (96) 


The expressions for H, and H, are now substituted into eq. 1, 


4 r 
along with eq. 6 to yield 
2 2 
V. ‘Rs 
- M1 2 A 2 





This is easily put into the dimensionless form used previously or 


De 
aes), = G2) [oY + (1 = 42)E mm 4)" (98) 


kq. 97 can now be substituted into eq. 88 and ® can be solved 
for straightforwardly. Iteration is necessary, however, since all loss 
ts depend on $, 
4.0 Loss Coefficient in Rotor 
Vavra (3) gives a relation for the losses in axial turbine 
bladings as 


Bronco Ze 2aee) fy lee ln ig | (99) 


104 
where AB is the flow deflection angle in degrees. The quantity Y is 
the velocity coefficient defined as the ratio of the actual and 


theoretical discharge velocities, which for the rotor is 


= (100) 
Ve eee 400 
Ro Wom 
anid 
AB, = Bs - B, (101) 
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The theoretical case occurs for Y, equal zero which in 


R 


pecordance Loe ed) somos vec treme Substituting eq. 100 into 


1 Bothy 
eq. 53 yields 





2 
ai | aneehs 22s peed 
Ouel Pp 9 2 ye 


The actual rotor blading losses using eq. 102 can now be written as 


Pp W,° p Wo° 
fae cae) = (+ 
2 ¥ 
Ww, /2 


In accordance with Horlock (Ref. 4) the loss coefficients depend 


on Reynolds numbers R,, and aspect ratio A, of ea. 60. Hence let Y 


R R 
% 
‘of eq. 103 now denoted as YR be the loss coefficient across the rotor 
= 
blading for Ap = 4 and Roe = 10°. The latter term is deIined as 
bes | 


yp = | (104) 


where DR is the hydraulic diameter of the discharge area of the flow 


channel between adjacent blades. It, is the rotor blade svacings at. 


o8 and Bn, is the opening at the discharge as shown in 


Figure 3, an approximation of the discharge opening is given by 


the radius R 


a, ~A, cos B, (105) 


then 


a: (flow area) _ 4 a, De (406) 
hR wetted perimeter 2 (a, = Dp) o 
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substituting for ap and rearranging yields 


ee (107) 
ve cos 8 
ao 2 


The blade spacing 4, can be found with respect to the radial 


extension Re - Ry by using the blade loading criterion af Zweifel (Ref. 5). 


For axial bladings this approximation is 


a 
— 2 (108) 
2 1 (Tan Bo - Tan 83) cos Bo 


Substituting this equation into the expression for the hydraulic 


diameter yields 





o( 8) 
D.. \z 
= r ba/® Tanb, = Jando j)coss (109) 
ea? EE B4p/COSP 5 


0,425(1-R, /R.,) 


Placing the above equation into the expression for the Reynolds number 
and recalling that 


OR, } 
Wo ~ Cos Bo (4) 


Sie Obtains 


b 
R 2 
2g) $ w Re, 


(b,/R,) (Lan8., - TanB,.) cos‘ B,, 
| cos, + | Vv 
: 0.425 Gaeta / 25 





Rup = 


(110) 


The loss coefficient is dependent on the aspect ratio and 


Reynolds number through the following relations after Horlock (Ref. 4). 


Yo eee (111 
R “2 ) 
! 1 AG 

Yp = ye (0915 ae ~ 1 (112) 


and 
10°\" ©? 


tp =i, _) (113) 


q 
where the coefficient YR holds for the actual aspect ration only at a 
Reynolds number equal to 10°, By combining the above equations the 


generalized loss coefficient is simply 


el 


= Dx ee) 
A a CaOr \ 
Yp = {op (.975 + PE) - ap GE) 
r= 2 \ eo aS 
These relations allow the actual loss coefficient of the energy absorber 
to be determined under varied conditions. 
Obviously the Reynolds number is not totally specified by 
geometric parameters, but depends on the flow coefficient ®. In 


iterating for ¢, the process must be assumed to begin at Y, , or a 


R 
Reynolds number of 10°, and subsequently correct the loss coefficient 
as improved values of $ and Rup are calculated. This procedure must 
also be employed for the velocity coefficient Y¥ which is dependent on 
@ through a dependence on the actual inlet angle Bae The iteration 


must be continued until convergence is achieved for all values 


dependent on 9. 
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4.1 Loss Coeiitelent igo catror 
This section follows closely the development in Section 3.0. 


Hag. 99 1S) aeain vwsed. with 


a 
- = (115) 
4Th 
and 
A = Of = iO 6 
Be Commie (116) 
For a frictionless process in the stator, the total pressure 
' t . 
(Pig) on at stator exit would equal total useful inlet pressure Pog 
or 
7 5 
(Pp_) P Ey ey) 
a (117) 


Brubeet puting eq. 115 into eq. 11/ yields 


e Pp V ‘ 
BD es ne (118) 
ay 2 
Ss 
For an actual flow 
2 
iF Pp V 
H es AE oils (119) 
a 
which yields through eq. 84 
1 es . 
‘ee ee (120) 
Ss 


This was as expected from Section 3.0. The loss coefficient can now be 


weroten after eqs. 112, 115, and 114 as 


a . | 7 
Y= 2 (.975 + — ; =f (121) 
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= 
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Y= oe (oe (22) 


“Ns 


and finally 


rg = 1 (975 + 928) - 4} ey” (123) 


where the aspect ratio is defined by eq. 90 as 


The Reynolds number for the stator is 


+ | (124) 


v 
where for a. ™ £7 Cosa » the hydraulic diameter is 


2b 


a 


Ne 
hs bY, 


1+——— 
Cos &@ 


4 


The blade spacing according to Zweifel (Ref. 5) is 


a | 0.425 
a? A (| Tan a, - Tan a | cosa 


4 5B. 4 


(126) 
and from eqs. 25 and 26 
y,=:o9R = = Gen 


The Reynolds number can now be expressed using the above four equations 


as 


qa 


‘Dp 
2 =) “Be @ w Ry, A 
pe 
Rus ~ bo D - 
we gee 5). Tan a, - Tan a, |) cos a, (128) 
R 
4 
0.425 (R,/R,) (1 - =) 


A series of iterations is again needed to insure proper 
convergence of all parameters, 
4.0 Pressure Loss Coefficients in Upper and Lower Passages 
In accordance with eq. 76 the pressure losses in the upper 


passage are expressed by 


(Ps), ~ Ps. Ys 
on (129) 


2 
(p.)m - P V 
le ey 1 M1 
ge ae oe (0) 


With Apb being the pressure drop in the bend, the loss 
coefficients of these bends defined as 
negro (131) 


are shown in Figure 15 as functions of the Reynolds number Bigs where 


similar to eqs. 104 and 124 





Ry = 2 (132) 


Vavra (Ref. 2) has developed the set of curves shown in Figure 15. 


These curves are based on the data of Sprenser (Ref. 6). Interpolation 


We 








was to obtain values for ratios of o/b different than those found 
in Sprenger. 


The hydraulic diameter D, is again 


_ 
D = 4 flow area (133) 


i wetted perimeter 


For a ring channel with flow area of 2m R b the wetted perimeter is 


2(2mR), ox 


2 iD 
th cae 4 ee (134) 

femece at the rotor inlet there is 
D, = 2b, (135) 


aagea toe stator inlet 


Dy, = 2b, (136) 


These above two values are introduced into eq. 132, with V 


equal to the corresponding value of Vins that iS, Vi before the rovor 
and Vung etienge: Gwicie: Sunsanena, VUMere) acesihUlie aus 
‘ V. 2b 
M4 S 
aaa Gion 
Vee © 
R, _ M1 R (138) 
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The Reynolds number is then used to obtain the loss coefficients § of 


Figure 5% 


j 


The loss coefficients € of Figure 15 have been obtainea from the 


t 


measured pressure drops by subratcting the pressure drops that would 


(2 





OCCuUC saminda straight duct at the corresponding Reynolds number and 
surface roughness. 

There are several conclusions which can be drawn from Figure 15 
and from additional data in Sprenger (Ref. 6). In general, the loss 
coefficients for bends made up of two concentric radii increase with 
decreasing *i/b ratios. The curves versus R, that exhibit a 
minimum indicate the flow separations occur at the inner radius, and 
should be avoided. Finally, Sprenger (Ref. 6) shows that two bends 
with Ti/b = 0.5 arranged in series to form a 180° bend have a loss 
coefficient of 1.66€ at Ry = 10° where € is the loss coefficient 
of one of the 99° bends. This data will be a first approximation in 
choosing the loss coefficient for the 180° bends in the energy 
absorber. 


4.1 Estimate of the Pressure Loss Coefficients in the Unver 





and Lower Passage 


The losses in the upper and lower passages will be assumed 
due to two components, ere bendinewand friepion. The Logeses due 
to duct bending were discussed in the previous section. An arbitrary 
assumption is to take these loss coefficients emtial! Tome, 0 Mancina) 
for the upper and lower passage respectively. The values Or ee vace 
given in Figure 15 for various value of my Asy it is recommended that 
ia/> be chosen between 0.5 and 1.0. Any bend with a value of ri/b 
larger than 1.0 will be assumed to have the same losses as a bend with 
*i/b equal 1.0. 


To the above losses must be added the frictional losses in a 


straight duct of hydraulic diameter D, and average length L, where 


h 
L; + L, 
L= oe (139) 
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and L. is the distance along the inner contour and L, that along the 
outer contour in the meridional plane from the trailing edges of one 
blade row to the leading edges of the other row. (See Figure 2) 

In accordance with eqs. 132 and 136 this general frictional 


pressure loss is 


bye (140) 


oe 


The friction factor }, is a function of both Reynolds number 


i 


and surface roughness. The value of hp is obtained from the curve of 


Figure 16, also from Vavra (Ref. 2), for a relative surface roughness 


of 
K 
i 1 
= (141) 
th 10 


where K. is the peak-to-valley roughness of the surface. Hence, for 


this case 


or > = 200 bd uw - in, (142) 
= iio 


If it is assumed that the rms roughness be one-half Kw the surface 
finish must be better than 100 b mnicro-inches. 

The applicability of this data to the three dimensional flow 
channels of the energy absorbed is obviously questionable. The losses 
G@erined below must , at best, be considered as first estimates only. 


The loss coefficient C. of eq. 129 will be taken as 


i 
as au 
C= 2.0€ + ab he (143) 
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C,=2.5§ + Bp” he (144) 


where € is obtained from Figure 15 and j, from Figure 16 for 


ea R er 
s M3 a Ky0 R Ss 
= ———— = (¢HwR = = =) — (145) 
Ryu . ( oe 5) 9 
and 
2b. V R 2b 
R ‘mM eae 
ee (146) 
Rye v ( 2k z) v 


In eq. 4 the length L, is given by 


, (lL), + (a), 
Ds aioe aes 


where (i). and Ca are the lengths ot the inner and outer 
Meridional contours, respectively, of the walls of the upper passage 


between stations 2 and 3. Similarly in eq. 133 


Cee a 
. . va ko 
a 5 | (147) 
where (L,). and (Ly), are the. counterparts of ee and (Li), for the 
lower passage between stations 4 and 1. 
5.O Losses of Angular Momentum in Upper and Lower Passages 
The frictional moments along the upper and lower passages 


were introduced by eqs. 18 and 22 respectively as 


Ma, =, ™Ro Vip | (148) 

and 
eS cs ho 
ee) ie al (149) 
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No data is presently available for the losses of angular 
momentum in passages of the type needed between the rotor and stator 
of the present absorber. A very crude approximation of these moments 
is obtained, however, by assuming that the shear stresses Me on the 
walls are merely those acting on a flat plate with a peripheral velocity 
Ve at R being the velocity outside the boundary layer. This over- 
simplification ignores the effect of the meridional velocity components 
which may cause large increases in the frictional moments. 

For fully turbulent boundary layers along rough walls the 


shear stress is 
= (ee 
Ve V (150) 


where Cr is assumed to be invarient with the Reynolds number of the 
Tiow. On the element of wall area en Rd L of Figure 2 there acts a 


frictional moment 


GM= = 2m hdl, jo 
ie O 


which becomes when using eq. 3 
(R al aids (151) 


It is assumed that the passages have been designed such that the flow 
is not separated from the walls, and all stresses external to the 
boundary layer are an order of magnitude smaller than stresses along 
the walls, and therefore negligible. This results in the product (RV, ) 


being constant outside of the boundary layers and equal to R the 


oY yo? 
value at the rotor discharge, for the upper passage. Integrating 


eq, 151 from station 2 to station 3 yields 


i 


(3) 


Mews : hia Tn ee) dL, (152) 
Me, = 70 Cp (Ry Vio) [(,), +), | (153) 


where the inner wall is of length (Ls and the outer wall is of 
length (Lo Hence by eq. 148 


pep Ry Yo | yy + (Lo | 


My (15!) 


mn 


Substitution of the mass flow equation and velocity relations into 
eq. 154 gives 


1+ 8 Tan B, l(t), = (1), | A 
Mo Op pee BS (155) 


Pl 


: Bo Oneaae 
from eq. 19 


1+ § Tan B, [(h,); +), | ye + (ZL ) al . 
Tn = + - %¢ a 2 be “2b Ky ee 


For the lower passage a Similar procedure and derivation leads 
to 
r 
+ 
Cp Tan a, | (Ly), (Ly), | 
My op a A SS (157) 


‘Bh S 


and 


By ee Cp aon Oy | (oy) ij (b,),| (158 ) 
2 Kay D 


80 





In eqs. 155 through 158, the absolute value of the flow angles 
must be introduced since the loss is a function of absolute deflection. 

For fully turbulent boundary layers, as in the present case, 
une local skin friction coefficients Ch depend on the surface roughness 


ration K/x; where K. is the peak to average value roughness and x is 


length along the plate. 


R 


1O 


I 


.OO4 for K_/x 
and 


1073 


R 


Cc 


, .006 for K_/* 


Hince there will be flow irregularities and possible flow -separations, 
po 8. = .016 was chosen 


since it gives best agreement with M-21 test data. 


it is advisable to use a larger value for c 
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APPENDIX B 


DETATLED ANALYSIS OF THE ENERGY ABSORBER WITH 
THE FLUID CIRCULATING IN THE REVERSE DIRECTION, MODES 3 AND 4 


A separate set of equations must be derived for the reverse flow 
modes. The reason, mentioned earlier, is that the relative flow exit 
angles which are assumed constant are now different. 

This development follows very closely to the previous analysis. 
Therefore, the intermediate steps of this derivation are omitted. 


1. Assumptions 


Bie By =—e@Mcuaae 
OQ = constant 
The remaining assumptions are identical to the normal flow case. 
2. Analysis 


The geometric relations employed are obtainable from Figures 6 


t 
and ‘. A new flow coefficient $ is conveniently defined as 


(B® i a 


The torque exerted on the fluid by the rotor when acting in the 


Gareccion of rotation is 
eo (cee eee ee) (2) 


where the mass flow equations are identical to those for the normal 


Case. 
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The frictional moments acting from stations 1 to 4 and stations 


3 to 2 are respectively. 


ee ae nl 


= 
I! 


Tu hy = R Oe 
where 


m R, V0 = By ae ~ Mey 


The relative flow inlet angles are 


aes ye 
Tan B, = "hen 5 is Tan a - @) x Ve 


S 


and 
(1 + 6 Tan 6, ) 


, b 
Ka Sk 


es 


Tan Sh 


The torque is now reduced to 
b 
: <) ; 
See Oa 5 Beianer 


where 


“pl 


b 
2 R 

Tt. = + 6 Tan B. -=— — Tan a 

R ; ab: 3 b. Wea a 
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(3) 


(4) 


(5) 


(6) 


(7) 


(8) 





Proceeding to the rotor, one can write, 


=, 
| 


_e Cos 1, 


R 7 Pop ~ Po 


_ oer 
ae 
pea 


and finally 











= Je ae! 
Ze 2 ie = I 
w R i 
an 
b 
: al R a 
a © {Tan Bi. - (= —— Tan @ ) cos ale } 
1 K3 bd. "Wa 3 R (9) 
reo < By) S. Pee ee 5 tR 
-° t D =) i( ) + (= Fn, Ten o,) | 4 
2 Kpo ‘ha 3 2cos By 
where 
ce i ee 
ee: 
Ro By 
for this case 
ANH . 
lA y = 33 (10) 
Ww RS 
In the lower passage 
Dp =e 
oe 
Oo Mt 
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and the energy drop is 


ee ee a 
2 -_— 7 


wR 4 
+ 8! 1( - we) Tan Bt (12) 


Oey ie Dae 2 (CRN ge ane : 
0? A EB ty. Ee Ge) na’, 


Net tor the stares 


II 


Vv 





4 vA 28 i 
1. Sep A 
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and the energy drop across the stator is 
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In the upper passage 


_ oma 7 Fe 


H, - H (eae. 2 Z b a2. 
—— = = E) Ca 7 cite Le zs aa o ) j 


w* R, e Kee ES 
The loss coefficients are now given as 
ae 
Ss 
ae 
R 
A., = ei o 
u eat ae = 
Ly 
A, =c, Tan 8, 
£ fi 1 bp Kay 


The Reynolds numbers ot obtain & and A, graphically are 


2b R, - ie 
a Sar G i i.) 


ug 


uae | pea 


(14) 


(i>) 


(16) 


(17) 


(18) 


The expressions for Y, Los and. x. remain the same as in the previous 


analysis, but the Reynolds numbers are now 
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fila 
4 Vv 
F ‘lee 
(d)/R, | Pan B, - Tan a cos’ B, 
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ase 


cosB, 4 


(19) 
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i w R 
2) 2) + 


D a 
Ss 2 
x [te OB -~ Tan a, | cos oe 


Rvs = 


cosa. + 


d 
Ry, ® 


0.425 (e) Cr - 1.0) 


The equation which is finally solved to obtain 6 is 


= H 
A He A Hy + A H. fu Hy 


eee oar (H, ~ H) = (H, ~ H) - CH, ~ H,) = 0 (20) 


The equation is solved by the program. An iteration technique must be 


{ 
carried out until all parameters dependent on ® converge. 
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APPENDIX C 


PROGRAM NORMAL 


Program NORMAL is essentially a modification of the program 
developed by NAEC. 

Program NORMAL solves the equations governing normal circulation, 
Modes 1 and 2, for a given design. The solution is based on an 
iteration technique which is repeated until the flow coefficient 
@ converges to within a specific accuracy. A starting value of 4 
is picked as are initial values of the various Reynolds numbers which 
are also %¢ dependent. The initial iteration step occurs with all 
losses based on standard initial Reynolds numbers. This subsequent 
woe Of 2 1S then used te improve the estimates of the Reynolds 
numbers and consequential losses, The technique is carried out 
until ¢ converges to the proper value. 

If the design parameters and angular speed are input data, the 
program will compute the flow coefficient, torque, velocities, losses, 
energy changes over each component, and the dimensionless torque, 
power, speed, and energy. The flow chart of Table IV depicts this 
process. 

The rated of data input is clearly indicated in the program 
itself. The program is felt to be self-explanatory in that comment 
cards have been extensively used. The Fortran symbols employed are 
given in the list of symbols. 

Subroutines have been used liberally. kach subroutine accomplishes 


one individual task such as determining a loss coefficient, or 


88 





enthalpy change across the stator, etc. The subroutines Fig. 15 and 
Fig. 16 are closest point approximations of the curves depicted in 
Figures 15 and 16 of the analysis. Thses points are given by Block 


data which is called by each subroutine. 
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APPENDIX D 


PROGRAM REVERSE 


Program REVERSE solves the equations for reverse circulation, 
Modes 3 and 4. The program is very similar to program NORMAL, and 
the flow chart of Table IV is the same. Further discussion would be 


mere repetition of Appendix C. 
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